Calculation of transition probabilities and ac Stark shifts in two-photon laser 

transitions of antiprotonic helium 



Masaki Hori 

Max-Planck-Institut fur Quantenoptik, Hans-Kopfermann-Strasse 1, 85748 Garching, Germany and 
Department of Physics, University of Tokyo, Hongo, Bunkyo-ku, Tokyo 113-0033, Japan 

Vladimir I. Korobov 
Joint Institute for Nuclear Research 141980, Dubna, Russia 

Numerical ab initio variational calculations of the transition probabilities and ac Stark shifts 
in two-photon transitions of antiprotonic helium atoms driven by two counter-propagating laser 
beams are presented. We found that sub-Doppler spectroscopy is in principle possible by exciting 
transitions of the type (n, L)—>(n — 2, L — 2) between antiprotonic states of principal and angular 
momentum quantum numbers n ~ L — 1 ~ 35, first by using highly monochromatic, nanosecond 
laser beams of intensities 10 4 — 10 5 W/cm 2 , and then by tuning the virtual intermediate state 
close (e.g., within 10-20 GHz) to the real state (n — 1, L — 1) to enhance the nonlinear transition 
probability. We expect that ac Stark shifts of a few MHz or more will become an important source 
of systematic error at fractional precisions of better than a few parts in 10 9 . These shifts can in 
principle be minimized and even canceled by selecting an optimum combination of laser intensities 
and frequencies. We simulated the resonance profiles of some two-photon transitions in the regions 
n — 30-40 of the p 4 He + and p 3 He + isotopes to find the best conditions that would allow this. 

PACS numbers: 36.10.-k, 31.15.A-, 32.70.Jz 



I. INTRODUCTION 



The transition frequencies u exp 



of antiprotonic helium 
atoms [1H5J] (pHe + = p~ + e~ + Hc 2+ ) have recently 
been measured by single-photon laser spectroscopy to a 
fractional precision of ~ 1 part in 10 8 |6|-|8|- By com- 
aring these results with three-body QED calculations 
the antiproton-to-electron mass ratio has been 
determined as 1836.152674(5) @, 0]. To further in- 
crease the experimental precision on is cxp , we have pro- 
posed future experiments [l5j of sub-Doppler two-photon 
spectroscopy of pHe + by irradiating the atom with two 
counter-propagating laser beams [16| . Dynamic (ac) 
Stark effects are expected to become one of the important 
sources of systematic error in these future experiments, 
as is the case with other high-precision laser spectroscopy 
measurements of atomic hydroge n I17H221 and antihydro- 
gen dHim, molecular hydrogen |25ll26|. helium f27H30lj . 
and muonium [3lH33j |. In this paper we calculate the 
transition probability and ac Stark shift involved in these 
two-photon transitions using precise three-body wave- 
functions of pHe + . 

ThepHe+ atoms @j can be easily synthesized b y si mply 
allowing an antiproton beam to slow down [34j-l39| and 
come to rest in a helium target. Some of the antiprotons 
are captured (40l - l49| into Rydberg pHe + states with large 
principal (n ~ 38) and angular momentum [L ~ n — 1) 
quantum numbers that have microsecond-scale lifetimes 
against antiproton annihilation in the helium nucleus. 
The longevity is due to the ground-state electron in pHe + 
which protects the antiproton during collisions with other 
helium atoms [H(3-[52j . All laser spectroscop y ex periments 
Q reported so far have used pulsed lasers [53} to induce 



single-photon transitions of antiprotons occupying these 
metastable states, to short-lived states with nanosecond- 
scale lifetimes against Auger emission of the electron [54]- 
[58j . A Rydberg pHe 2+ ion [59rl62T | then remained after 
Auger decay, which was rapidly destroyed by collisional 
Stark effects. The resulting resonance profiles of pHe + 
had Doppler widths Alj£>/2tt ~ 0.3 — 1.2 GHz corre- 
sponding to the thermal motion of pHc + in the experi- 
mental target at T = 10 K. This broadening limited the 
experimental precision on ^ cxp . 

The first-order Doppler broadening can in principle be 
reduced (Fig.[T]) [15j by irradiating pRe + atoms with two 
counter-propagating laser beams of angular frequencies 
oji and u>2 and inducing, e.g., the two-photon transition 
(n, L) — > (n — 2,L — 2). This results in a reduction of 
Atop by a factor |(o;i — u)2)/(uj\ -\-uj2)\- Among a number 
of possible two-photon transitions, a particularly strong 
signal is expected for (n, L) = (36, 34) -> (34, 32UFig. [TJ, 
as this involves a large antiproton population [40, [y| m 
the resonance parent state (36,34). Whereas the states 
(36, 34) and (35, 33) are metastable with 1-^s-scale life- 
times, the resonance daughter state (34, 32) is Auger- 
dominated with the lifetime r ~ 4 ns [y, [l(| 0, HH that 
corresponds to the natural width of a spectral line of ~ 30 
MHz. 

The probability of inducing the transition can be en- 
hanced [l5l . [63Tj63 | by tuning wi and 0J2 so that the virtual 
intermediate state of the two-photon transition lies close 
(e.g. \Auj d /2ir\ < 10-20 GHz) to the real state (35,33) 
such that, 



Wi(SI) = 47ri? 00 C [#(35,33) - #(34,32) (a-u.)] + Au d , 



w 2 (SI) = AnRooC [E { 



(36,34) — #(35,33) 



,.U.)] - Acjrf. 



(1) 
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Here R^c = 3.289842 x 10 15 Hz denotes the Rydberg 
constant, and the binding energy of the p 4 Hc + 

state (n, L) in atomic units. At experimental conditions 
wherein this offset is much larger than the Doppler width 
(|AW(j| 3> Awe), the two-photon transition is expected 
to directly transfer the antiprotons populating the par- 
ent state (36,34) to the daughter state (34,32), whereas 
the population in the intermediate state (35, 33) will be 
unaffected. 

This paper is organized in the following way. Some 
details of the numerical methods are described in Sec- 
tion [Til The transition amplitude of the two-photon res- 
onance (36, 34)— >-(34, 32) of p 4 He, and the polarizabilities 
of the parent and daughter states (36, 34) and (34, 32) at 
various offsets Aojd are estimated in Sections IIII Al and 
IIIIBI We next calculate the "background" polarizabil- 
ity due to the contributions of pHc + states other than 
the resonant intermediate state (35,33) ( Section IIII CI) . 
Based on these results, the ac Stark shift and broaden- 
ing are semi-analytically estimated in Section IIII Dl We 
next discuss the hypcrhnc structure in two-photon tran- 
sitions of the p 4 He + and J5 3 He + isotopes (Section IIII Ep 
[76T - l8ll ]. We numerically simulate the profile of several 
two-photon resonances (Section IIII F[) before concluding 
the paper. Atomic units are used to evaluate the state po- 
larizabilities and transition amplitudes (70l - [75l |. whereas 
International System of Units (SI) are used for the tran- 
sition frequencies, rates, and laser intensities relevant for 
future spectroscopy experiments. 



II. DETAILS OF THE CALCULATION 

For simplicity we take the linearly polarized laser field 
aligned along the z-axis, such that the perturbation 
Hamiltonian in a laser field of frequency uj and ampli- 
tude F has the form. 



H' = -e¥(t)d, F(t) = e z Fcos(wt). 



(2) 



Here d = Yl a =i ^aR-a is the electric dipole moment op- 
erator. The second order correction E^ to the unper- 
turbed eigenenergy Eq of a pHe + state vector |0) may 
then be expressed as, 



£( 2 ) = -^af( W ,M 2 ).F 2 , 



(3) 



where oy* (w) is a tensor of the dynamic dipole polariz- 
ability: 



E 



(0|<%)<g|^|0) , {0\#\q)(q\d?\Q) 



Eq — E q + to Eq — E q 



(4) 



The energy of a pHe state vector \q) is denoted by E q , 
and the summation of q is over all states which are ac- 
cessible via single-photon transition from the resonance 
parent state |0). 



The tensor a^(uj) may be rewritten in terms of the 
irreducible scalar and tensor polarizability operators, 



a i 



a + + a + O-l , 



(5) 



a+L(2L-l) o (2L-l) 



3(Z+l)(2Z+3) 3(i+l) 



where the angular momentum operator is denoted by L. 
The coefficients a + , ao, and a_ are defined as follows, 



a+ 



ao 



(£; -f;,)|(OL||d||g(L + l))| 



2L 
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{ E -E q )\(OL\\d\\qL)\ 

(E - E q y - ^ ' 
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(E - E q f - W 2 



(6) 



Here a+ and a_ represent the contributions from antipro- 
ton transitions to states of normal parity which change 
the orbital angular momentum quantum number L of the 
antiproton by 1 or — 1 . The contribution ao involves tran- 
sitions to pHe + states of anomalous parity wherein the 
L-value is unchanged and the ls-electron is excited to, 
e.g., the 2p state. 

For our analysis it is convenient to define " background" 
polarizabilities using the above equations, wherein the 
dominant contribution from the intermediate state of the 
two-photon transition is subtracted. For example, the 
corresponding scalar and tensor background polarizabil- 
ities of state (n, L) — (36, 34) can be calculated as, 



at 



2(E -E l ) |(0L||d||i(L-l))| 
3(2£+l) 
2(E -Ei) 



{Eo - Ei) z 
|(0L||d|K(L- 



!))| 



(7) 



3(2L+1) (Eo-Et) 2 



Here E and Ei denote the energies of state (36,34) and 
the intermediate state (35,32). 

The transition matrix element kl,l-2,m of the pHe + 
two-photon transition (n, L, M)— >(n— 2, L— 2, M) induced 
by two linearly-polarized laser beams can be calculated 
using the Wigner 3j-symbols as, 



L-l 1 L 

M -M 



L-l 1 L-2 
M —M 



(8) 



(l(X-2)||d|| g (L-l))( g (L-l)||d||0L) 
Eo — E q — loi 

| (l(£-2)||d|| g (£-l))( g (X-l)||d||0L) 
Eo — E q — u>2 

wherein |1) denotes the state vector of the resonance 
daughter state. This is related to the two-photon Rabi 
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oscillation frequency (in atomic units) of this laser tran- 
sition via the equation, 



= ~ \k L ,L-2,M\ F 1 F 2- 



(9) 



The last term in Eq. ([8]) gives a small contribution and 
will be neglected. 

In order to calculate these quantities we must evaluate 
the reduced matrix elements for the dipole operator d 
and diagonalize the Hamiltonian. For this we use the 
variational exponential expansion described in Ref. 
The wave function for a state with a total orbital angular 
momentum L and of a total spatial parity n = (— 1) L is 
expanded as follows, 

*I A/ (R,r 1 )= y, y^UA^Gfj^n,^), 

h+h=L 
N 

Gfj 2 (i?,n,r 2 ) = ]T {^Re^-^ 1 ^^] ( 10 ) 

71=1 

+AJm[e- a " fl " f, " ri - 7 " r2 ]}. 

where the complex exponents a, /3, and 7 are generated 
in a pseudorandom way, R and ri are position vectors of 
an antiproton and an electron with respect to a helium 
nucleus, and ri the distance between the antiproton and 
electron. Further details may be found in Refs. [9l. Hoj. 

III. RESULTS 

A. Transition matrix element 

We first calculated the transition amplitude kl,l-2,m 
for the two-photon resonance (n, L) — (36, 34)— >-(34, 32) 
in |5 4 He at various offsets uid from the intermediate 
(35, 33) state, and estimated the laser intensities needed 
to drive this transition. In Fig. O sequences of single- 
photon transitions connecting the states (36,34) and 
(34, 32) are indicated by straight arrows, together with 
the corresponding dipole moment (0L||d||g(L — 1))|. The 
atomic units shown here can be converted to SI units, the 
corresponding spontaneous decay rates in s _1 may be ob- 
tained using the equation, 



7 (SI) 



e 2 a 2 4u, 3 |(Q£||d|| g (£-l))| 2 (a.u.) 
47re 3?ic 3 2L + 1 



(11) 



The Si-unit constants that appear in the above equation 
are, e: the elementary charge, ao'. the Bohr radius, eo: 
the dielectric constant of vacuum, Ti: the reduced Planck 
constant, and c the speed of light. The angular transition 
frequency between states \q) and |0) are denoted by tu q o. 

Two types of transitions (n, L) — > (n,L— 1) and 
(n, L)— >(n— 1, L— 1) have the largest amplitudes of ~ 1 
a.u., but the later kind which conserve the vibrational 
quantum number v = n — L — 1 and involve fluorescence 
photons of frequency Lo q [)/2iT ~ 10 15 Hz constitute the 



dominant channels of spontaneous decay. These tran- 
sitions are most favorable for laser spectroscopy. The 
transition frequencies, dipole moments, and decay rates 
of some single-photon resonances in p 4 He + and |5 3 Hc + of 
the type (n, L)— >{n— 1, L— 1) are shown in Table|U The 
dipole moments for the higher-lying infrared transitions 
involving states with n ~ 40 are relatively large (> 2 
a.u.), whereas for UV transitions in the n < 33 regions it 
is reduced to < 1 a.u. On the other hand, the radiative 
decay rates increase for lower-ro transitions, e.g. from 
4.7 x 10 5 s- 1 for (40,36)^(39,35), to 6.6 x 10 5 s" 1 for 
(32, 31)— ^(31,30) due to the c<j 3 -dependence. 

Using the single-photon dipole moments calculated 
above, we derive the two-photon transition amplitude 
k-l.l-2,m of the resonance (36, 34)— >-(34, 32) in p 4 Hc + for 
cases wherein the virtual intermediate state is offset over 
a large range between Aw(i/27r = —0.7 and 0.7 PHz from 
the state (35,33). The atom is excited by two linearly- 
polarized, counterpropagating laser beams. Fig.|3](a)-(c) 
show the amplitude k,l,l-2 averaged over all ~ 70 transi- 
tions between the magnetic substates which conserve the 
M-value. The |K£.L_2|-values are usually small, e.g. ~ 1 
a.u. for lasers of equal frequency (w\ — L02). This is an 
order of magnitude smaller than the amplitude ~ 7.8 a.u. 
[20j for the ls-2s two-photon transition of atomic hydro- 
gen excited by 243-nm laser light. Gigawatt-scale laser 
intensities would be needed to induce the antiprotonic 
transition within the microsecond-scale lifetime of pHc + . 
On the other hand, the transition probabilities can be 
strongly enhanced to > 10 3 a.u. by tuning the virtual 
intermediate state within ~ 20 GHz of the real states 
(n,L) = (34,33), (35,33), or (36,33). Eq. U indicates 
that the transition can then be induced using nanosec- 
ond laser pulses of electric field F ~ (1 — 2) x 10~ 6 a.u. 
According to the equation, 

/(SI) 



1 /eF(a.u.) 
2 £oC 



(12) 



this corresponds to a peak intensity of / ~ 10 4 -10 5 
W/cm 2 which is achievable using titanium sapphire lasers 
of narrow linewidth II 611 . 



B. Polarizabilities 

We next evaluate the polarizabilities of the parent 
and daughter states of the transition. In Figs. |4] (a) 
and (c), the scalar and tensor polarizability components, 
oe s (wi) (34,32), a* (^i) (34,32), of state (n,L) = (34,32) are 
shown. To simplify the calculation, we initially assume 
that the atom is irradiated by a single laser field of fre- 
quency loi (see Fig. [IJ corresponding to offsets between 
Aoj d /2n = -20 and 20 GHz from the state (35,33). 
A similar plot for the polarizabilities Q: s (oj2)(36,34) an d 
c*t(w2)(36,34) °f state (36, 34) irradiated by a laser field of 
u>2 are shown in Figs. [4] (b) and (d). 

As the respective lasers are offset from Auid/^iT = 
— 100 to -6 GHz, the scalar polarizabilities decrease from 
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~ —500 to ~ — 1 x 10 4 a.u., whereas the tensor polariz- 
abilities have opposite sign and increase from ~ 500 to 
- 1 x 10 4 a.u. (Table HTJ)- Polarizabilities of 1000-10000 
a.u. correspond to 50-500 Hz/(W/cm 2 ) in SI units ac- 
cording to the equation, 



a(a.u.) = a(SI) x 



(13) 



These graphs follow a reciprocal AwJ dependence and 
are approximately symmetric with respect to the origin. 
This simple behavior suggests that the ac Stark shift is 
primarily caused by the contribution from the interme- 
diate state (35, 33). 



C. Contribution of nonresonant states 

We next study the non-resonant or " background" con- 
tributions (Eq. [7} to the polarizabilities from all pHe + 
states other than the intermediate state (n, L) = (35, 33). 
This will allow us to estimate how far the measured ac- 
Stark shift would deviate from the predictions of a sim- 
ple three-level model which include only the resonance 
parent, daughter, and intermediate states. Figs. [5] (a) 
and (c) show the background scalar and tensor polariz- 
abilities P s {u2) (36,34) and p t (wi) (36,34) of state (36,34) 
when irradiated with a single laser field of frequency uj 2 
at offsets between Awd/27r = —20 and 20 GHz. They 
remained relatively constant at ~ —0.40 and 0.95 a.u. 
respectively (Table HT|). Figs. [5] (b) and (d) are the corre- 
sponding plots of /3 s (wi) (34,32) and Pt (^l) (34,32) for state 
(34, 32) irradiated with the wi-laser. They are simi- 
larly constant (—1.36 and —0.40 a.u.). All these back- 
ground polarizabilities are at least three orders of magni- 
tude smaller than the dominant contributions to a s and 
at arising from the intermediate state (35, 33) at offsets 
|Aw rf /27r| < 12 GHz. 

The case of two counterpropagating laser fields of an- 
gular frequencies uj\ and ui 2 , and amplitudes F± and F 2 
irradiating the atom simultaneously will next be consid- 
ered. The perturbation Hamiltonian for this can be ex- 
pressed as, 



W = -eF(t)d, 

F(t) = e z [Fi cos(w 1 t) + F 2 cos(w 2 t)} . 



(14) 



As shown in Ref. [17| , the interference effect between the 
two laser fields can be neglected in the case of u>\ / ui 2 . 
The contribution to the ac Stark shift of state (36, 34) 
from the wi-laser [which is far off-resonance with respect 
to the upper single-photon transition (36, 34)— >(35, 33)] 
is expressed by the scalar and tensor polarizabilities 
a a (u>i) (36,34) and at (u)i) (36,34)- The calculated values at 
offsets between AuJd/2TT = —20 and 20 GHz were re- 
spectively — 1.6 and —0.2 a.u. (Table |TTJ) . The corre- 
sponding values for the daughter state a s (uj 2 ) (34,32) and 
c*t(w2) (34,32) were also small (—2.1 and 0.3 a.u.). 



We conclude that the two-photon spectroscopy exper- 
iment depicted in Fig. Q] can be accurately simulated by 
a simple model involving three states interacting with 
two laser beams. Any non-resonant contribution from 
other pHe + states are at least three orders of magnitude 
smaller. 



D. ac Stark shifts of transition frequency 

The ac Stark shift in the transition frequency of the 
(36, 34)— s>(34, 32) resonance can be analytically estimated 



2vr 



2R a 



E, 



(2) 

(36,34,M) 



E 



(2) 

(34,32, M) 



(15) 



wherein the ac Stark shifts E, 



(2) 



r(2) 



and E (34,32,M) 111 



J (36,34,M) 

the parent and daughter states of magnetic substate M 
induced by the two linearly-polarized laser fields can be 
approximated as, 



E, 



(2) 

(36,34,M) 
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~2 
3M 2 



«s( w 2)(36,34) 
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2278 



a t {u 2 ) 



(36,34) 



Fo 
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E. 
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(34,32,M) 
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3Af 2 



tts(^l) (34,32) 
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2016 



£*t(u>i) 



(34,32) 
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In Fig. [6] (a), the ac Stark shift Aw ac /27r for M-values 
0, ±16, and ±32, and two laser offsets Au)d/2n = —12 
and 12 GHz obtained from the above equations are plot- 
ted as a function of the intensity ratio I 2 /Ii = F 2 /Fi 
between the two laser beams. Here 1\ fixed at ~ 5 x 10 4 
W/cm 2 while I 2 is scanned between (3 — 5) x 10 4 W/cm 2 . 
We find a positive shift (Acj ac > 0) at two combina- 
tions of laser offsets and intensities (Aw^ > 0, I\ ~3> I 2 ) 
and (Awd < 0, I\ <C Is)- Conversely the shift is nega- 
tive Aw ac < at (Au} d < 0, I\ > I2) and (Auj d > 0, 
I\ *C I 2 ). In addition to the ac Stark shift, the tensor 
polarizabilities cause the resonance line to split depend- 
ing on the M-value of the involved states. At conditions 
of I 2 /Ii < 0.6 or > 1 the ac Stark shift and splitting 
can reach values of more than 5-10 MHz. At smaller off- 
sets Au}d/2iT — ±6 GHz, the ac Stark shift and splitting 
become twice as large [Fig. [B] (b)]. 

The ac Stark shift arising from a s and the splitting due 
to at can be minimized by adjusting the laser intensities 
to the values indicated by filled circles in Figs. |6] (a)-(b), 



h 



a«(wi)( 



34,32) 



h a s (u 2 ) (36,34) 



0.76. 



(17) 



An important point is that when the sign of the laser 
detuning is reversed, e.g., from Aujd/2i: = —12 to 12 



5 



GHz, the resulting ac Stark shift also reverses sign but 
its magnitude remains constant, 

Aw ac (wi,o;2) = -Acj ac (wi -2Acj ( j,W2 + 2Aw ( i). (18) 

This means that the residual ac Stark shift can be can- 
celed by comparing the two-photon transition frequencies 
measured at laser offsets of opposite sign but the same 
absolute value, i.e., Aujj, and -Aw,;. 

Fig. [7] shows the ac Stark shifts and transition ampli- 
tudes kl,l-2,m of all magnetic sublines between M — 
—32 and 32 of this two-photon resonance at laser off- 
sets Aoj d /2n = 12 GHz (a)-(c) and -12 GHz (d)-(f). 
The profiles were calculated at three ratios of the laser 
intensities J2//1 = 0.6, 0.78, and 1.0 with 1\ fixed at 
5 x 10 4 W/cm 2 . The ac Stark effect causes a triangular 
profile with the M — transitions being strongest and 
shifting the most. In actual experiments involving pulsed 
laser beams, these shifts and splittings will smear out due 
to nanosecond-scale changes in the light field intensities 
h{t) and I 2 {t). 



E. Hyperflne structure 

We next study the hyperflne lines that appear in the 
two-photon resonance profile. The hyperfme substates 
of the parent, intermediate, and daughter states (n,L), 
(n — 1, L — 1), and (n — 2, L — 2) in |5 4 He + are shown 
schematically in Fig. [5] (a) . Due to the dominant inter- 
action between the electron spin S e and the antiproton 
orbital angular momentum L, a pair of fine structure 
sublevels of intermediate angular momentum quantum 
number F = L ± 1/2 and a splitting 10-15 GHz arise. 
The interactions involving the antiproton spin S-p- cause 
each fine structure sublevel to split by a few hundred 
MHz into pairs of hyperfme sublevels of total angular 
momentum quantum number J = F ± 1/2. In Fig. [5] 
(a), the spin orientations of the electron and antiproton 
(S e ,Sp) are indicated for the four hyperfme sublevels. 
For example the energetically highest-lying component 
consists of a spin-down electron and spin-up antiproton, 
i.e. (S e ,Sp) = (it)- 

In the J5 3 He + case [Fig.|S](b)], the electronic fine struc- 
ture sublevels of F — L ± 1/2 are similarly split by ~ 10 
GHz. Each fine structure sublevel is then split into pairs 
of 3 He hyperfme sublevels of intermediate angular mo- 
mentum K = F ± i arising from the interactions involv- 
ing the nuclear spin S/j. The antiproton spin gives rise 
to eight hyperfme sublevels of total angular momentum 
J = K±7j. The spin orientations of the three constituent 
particles (S e , Sh, Sp) are indicated for each substate in 
Fig. Kb). 

In Fig. |S] (a)-(b), the four and eight strongest two- 
photon transitions between the hyperfme sublevels of 
p 4 He + and p 3 He + are indicated by arrows. These tran- 
sitions pass through the virtual intermediate state with- 
out flipping the spin of any constituent particle. Many 



other transitions are possible, but they all involve spin- 
flip and so their amplitudes are suppressed by three or- 
ders of magnitude or more. 



F. Optical rate equations 

To simulate pHc + two-photon resonance profiles, we 
used the following nonlinear rate equations which de- 
scribe a three-level model, 



dp a 



dt 

dpbb 
dt 

dpcc 
dt 

Opab 



—lm(CliMPab), 

lm(fti M p ab ) - lm(n 2 MPbc), 

--fcPcc + lm(fl 2 MPbc), 



.^1M, \ , .^2M 



(19) 



—QjT = -idabPab + i—^—(Paa ~ Pbb) + i—^—Pac, 

dpbc . , , . &>2M / \ . 

= -idbcPbc + l—^—(Pbb - Pec) - l—j—pac, 

dpab ., , .^2M MlM 

—Qj~ = -idacPac + i—^-Paa ~ i—^~Pbc- 



Here the density matrix p aa , pbb, and p cc represent the 
antiproton populations in the parent, intermediate, and 
daughter states. The mixing between pairs of states in- 
duced by the lasers are denoted by p ab , p bc , and p ac , the 
Auger decay rate of the daughter state by j c . The three 
detunings that appear in Eq. Q1J] can be calculated using 
the equations, 



d a b — Eb — E a — 



dbc 
doc. 
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+ 7c 
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) Ui 




c - 



.la+lc 



(20) 



where v z denotes the velocity component of the atom 
in the direction of the wi-laser beam, E ai Eb, and E c 
the binding energy of the hyperflne states involved in the 
transition, and i a and 7b the radiative rates of the parent 
and intermediate states. The angular Rabi frequencies 
of single-photon transitions induced between the parent 
and intermediate states of magnetic quantum number M , 
and the intermediate and daughter states are respectively 
denoted by Q\m and f^M- 

Values of 57im/27t (in s _1 ) for the transition 
(0, L, F, K, J, M)^t(q, L—l, F' , K', J', M) in p 3 He + by 
linearly-polarized laser light of intensity l\ (in W/cm 2 ) 
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can be calculated, 



<ws.i.) = M^. 

V eoc h 



(21) 



(OLFKJM\d\q(L-l)F'K'J'M) (a.u 



where the transition matrix element (in atomic units) can 
be derived using the Wigner 3j- and 6j-symbols, 



; 0LFKJM\ d \q(L-l)F'K'J'M ] 



(oi||d|| ? (i-i))(^ I J M 



x^(2J+l)(2J' + l) 



A" J' i 
J A 1 



(22) 



*^{2K + l){2K' + l)\ 5 



Xv/(2F+1)(2F'+1) 



F L 1 



Eq. ([22]) only provides approximate values for the tran- 
sition matrix elements since L, F, and K are not exact 
quantum numbers of the three-body Hamiltonian. The 
results however agree well with exact transition ampli- 
tudes within 1% accuracy. 

We assume that the pHe + M-states are uniformly pop- 
ulated in the parent state and follow a Maxwellian ther- 
mal distribution with T~ 10 K. We integrate Eq. to 
simulate the antiprotons depopulated by the two laser 
beams into the daughter state producing the spectro- 
scopic signal. The laser pulses are assumed to have Gaus- 
sian temporal profiles of pulse length At = 100 ns. 

Fig. [9] (a)-(b) shows the efficiency e of the laser pulses 
depleting the population in the parent state of the two- 
photon transition (36, 34)— >(34, 32) of p 4 He + [i.e. e = 1 
if the laser induces all the antiprotons occupying state 
(n, L) to annihilate, and e — when no such anni- 
hilations occur]. The virtual intermediate state is off- 
set LOd/lir = —12 GHz away from state (35,33), so 
that uj\ + u>2 coincides with the hyperfine component 
(S e ,Sp) = (tt)-Ktt) 01 tne resonance line. As the laser 
intensity is increased between p = and 0.2 mJ/cm 2 , 
the e- value increases quadratically as expected for a two- 
photon process. It begins to saturate at p > 1 mJ/cm 2 
corresponding to e ~ 0.3. Monte Carlo simulations indi- 
cate that the two-photon resonance signal corresponding 
to e ~ 0.3 would be strong enough for clear detection 
against the background caused by spontaneously annihi- 
lating antiprotons [82J with a signal-to- noise ratio of > 5. 
Higher laser intensities would of course provide an even 
stronger signal, but power broadening effects then deteri- 
orate the spectral resolution to several hundred MHz and 
so this should be avoided for high-precision spectroscopy. 

The resonance profiles of the two-photon transitions 
(40,36) (38,34) of the v = 3 cascade, (38,35) 



(36.33) of v = 2, (36,34) (34,32) of v = 1, and 
(33,32) — > (31,30) of v — in p 4 He + at temperature 
T ~ 10 K are shown in Figs. Q1 (a)-(d). These res- 
onances have among the largest transition amplitudes, 
and the Auger decay rates of the daughter states are 
large ja = (2-5 — 4) x 10 8 s~ 4 which is a neces sary 
condition to obtain a strong annihilation signal [82j . 
The intensities of the two lasers are around p ~ 1 
mJ/cm 2 . The laser frequency uii is fixed to an offset cor- 
responding to Aujd/^TT — —12 GHz from the intermediate 
state, whereas u>2 was scanned between -0.9 and 0.9 GHz 
around the two-photon resonance defined by u>\ + ui2- In 
each simulated profile, the positions of the four hyperfine 
lines are indicated with arrows together with the cor- 
responding spin orientations (S e ,Sp). The ~ 200 MHz 
linewidth of these profiles are primarily caused by the 
large Auger width of the daughter states, and the resid- 
ual Doppler and power broadening. 

The resonance (40, 36)— >(38, 34) shows a two-peak 
structure [Fig. [10] (a)] with a frequency interval of ~ 1.1 
GHz which arises from the dominant spin-orbit interac- 
tion between S e and L. Each peak is a superposition 
of two hyperfine lines with a few tens of MHz spacing 
caused by a further interaction between the antiproton 
and electron spins. The asymmetric structure of the pro- 
file of Fig. [TO] (a) is due to the fact that the 25-MHz 
spacing between the hyperfine lines (S e , Sp) = (tt)-Ktt) 
and (tl)— Ktl) are small compared to the 75-MHz spac- 
ing between (It) - Kit) an d (ID-KII) The spacings be- 
tween the hyperfine lines becomes gradually smaller for 
lower-lying transitions involving states of smaller n- and 
w-values, e.g., 0.8 and 0.5 GHz for (38, 35)->(36, 33) and 

(36. 34) — K34, 32). The hyperfine lines can no longer be 
resolved for the lowest— n transition (33, 32)— >(31, 30) 
[Fig. rfj] (d)]. The low transition probability (Table [J) 
of this resonance causes the small depopulation efficiency 
£ < 0.1 seen here; laser intensities of p > 2 mJ/cm 2 would 
be needed to produce a sufficient experimental signal. 

We expect the two UV transitions (n, L) = 
(36, 34)^(34, 32) and (33, 32)->(31, 30) in p 4 He + to yield 
the highest signal-to-noise ratios in laser spectroscopy ex- 
periments. This is because the parent states (36, 34) and 
(33, 32) retain large antiproton populations for long pe- 
riods t = 3-10 ^s following pHc + formation [40] - By 
comparison, cascade processes rapidly deplete the pop- 
ulations in higher n > 37 states within 1-2 fis, and so 
the associated two-photon spectroscopy signals contain a 
large background due to the spontaneously annihilating 
pHe + atoms [toLElj. 

Higher experimental precisions on v cxp may be 
achieved by cooling the atoms to lower temperature 
and by inducing two-photon transitions between pairs of 
pHe + states with microsecond-scale lifetimes. Fig.[TU](c) 
shows the resonance (37, 35)— >(35, 33) of p 4 He + at tem- 
perature T — 1.5 K. Both parent and daughter states 
have lifetimes of r ~ 1 fis, and so its natural linewidth 
~ 200 kHz is two orders of magnitude smaller than in 
the other resonances Fig. [TU] (a)-(d) studied here. It is 
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unfortunately difficult to measure this transition exper- 
imentally, as the present detection method requires the 
daughter state to proceed rapidly to antiproton annihi- 
lation. 

The profiles of two p 3 He + resonances which are ex- 
pected to yield the highest signal to noise ratios 

(35.33) ^(33,31) and (33,32)^(31,30) at temperature 
T ~ 10 K are shown in Fig. [TTJ (a)-(b). The positions 
of the eight hyperfine lines and their spin configurations 
(S e , Sh, Sp) are indicated by arrows. Due to the large 
number of partially overlapping lines, it may be difficult 
to determine the ^ oxp -values for p 3 He + with a similar 
level of precision as in p 4 He + . The problem would be 
especially acute in the case of (33, 32)— >(31, 30) [Fig. [11] 
(b)] which contain 8 sublines within a relatively small 
0.6- GHz interval. 

We finally use these numerical simulations to deter- 
mine the ac Stark shift under realistic experimental con- 
ditions. Fig. [T2] (a) shows the profiles of the resonance 

(36.34) ^.(34,32) of p 4 He + at temperature T~10 K and 
laser offset Acj t ;/27r = — 12 GHz. They were calculated 
at two combinations of the laser intensities: I\ = 5 x 10 4 
W/cm 2 and I2 = 2.5 x 10 4 W/cm 2 (broken lines) and 
h = h = 5 x 10 4 W/cm 2 (solid lines). As I 2 /h is in- 
creased, the transition frequency shifts to larger values. 
In Fig. [TJ] (b), the ac Stark shifts Aw ac /27r determined 
from the simulated profiles of Fig. [T2] (a) at laser offsets 
Aud/2ir = — 12 GHz are plotted using filled circles. It 
increases linearly from —2 MHz at I2/I1 — 0.5, to 5 MHz 
at Iz/Ix = 1. A similar plot for offset Au>d/2n — 12 GHz 
is shown using filled squares. The two calculated sets of 
ac Stark shifts are of equal magnitude and opposite sign, 
the minimum occurring around ~ 0.65. 

IV. CONCLUSIONS 

We conclude that two-photon transitions in pHc + of 
the type (n, L)— >(n — 2, L — 2) can indeed be induced us- 



ing two counterpropagating nanosecond laser pulses of 
intensity ~ 1 mJ/cm 2 , for cases wherein the virtual in- 
termediate state is tuned within \AuJd/2ir\ = 10-20 GHz 
of the real state (n— 1, L—l). The spectral resolution of the 
measured resonances should increase by an order of mag- 
nitude or more compared to conventional single-photon 
spectroscopy. The ac Stark shifts at these experimental 
conditions can reach several MHz or more, but this can 
be minimized by carefully adjusting the relative intensi- 
ties of the two laser beams. Any remaining shift can be 
canceled by comparing the resonance profiles measured 
at positive and negative offsets ±Au)d of the virtual in- 
termediate state from the real state. In practice, this 
can be done by, e.g., using a frequency comb |83] to ac- 
curately control the frequencies oji and C02 of the coun- 
terpropagating laser beams. The UV two-photon transi- 
tions (36,34)^(34,32) and (33,32)^(31,30) in p 4 He + , 

and (35, 33)— J>(33, 31) in p 3 He + are expected to yield par- 
ticularly strong resonance signals that can be precisely 
measured. 
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TABLE I: Some single-photon transitions in p 4 He + and j^He : their vibrational quantum number v, transition frequency, and 
absolute value of the transition matrix element (i.e. dipole moment). 

(n, L) —¥ (n' , L') v Trans, freq. Dipole moment Rad. rate 









(THz) 


(a.u.) 


(10 5 s" 








p 4 He + states 






(40, 36) - 


->■ (39, 35) 


3 


444.8 


2.28 


4.71 


(39, 35) - 


->■ (38, 34) 


3 


501.9 


2.02 


5.44 


(38,35) - 


->■ (37, 34) 


2 


566.1 


1.82 


6.36 


(37,34)- 


■> (36, 33) 


2 


636.9 


1.58 


7.01 


(37,35) - 


->■ (36, 34) 


1 


638.6 


1.61 


7.15 


(36, 34) - 


->■ (35, 33) 


1 


717.5 


1.38 


7.63 


(35,33) - 


->■ (34, 32) 


1 


804.6 


1.16 


7.92 


(33,32) - 


->■ (32,31) 





1012.4 


0.79 


7.42 


(32,31) - 


->■ (31,30) 





1132.6 


0.62 


6.62 








p 3 He + states 






(39,35) - 


->■ (38, 34) 


3 


445.8 


2.30 


4.96 


(38, 34) - 


-> (37, 33) 


3 


505.2 


2.03 


5.77 


(37, 34) - 


->■ (36, 33) 


2 


572.0 


1.83 


6.81 


(36,33) - 


->■ (35,32) 


2 


646.2 


1.58 


7.55 


(35,33) - 


->■ (34, 32) 


1 


730.8 


1.37 


8.27 


(34, 32) - 


->■ (33,31) 


1 


822.8 


1.15 


8.59 


(33,32) - 


->■ (32,31) 





928.8 


0.95 


8.44 


(32,31) - 


->■ (31,30) 





1043.1 


0.77 


7.94 



TABLE II: Scalar and tensor polarizabilities of the p 4 He + states (34, 32) and (36, 34) at offset frequencies of the lasers between 
Au d /2ir = -100 and 100 GHz, see text. 



Polarizabilities (a.u.) 



Acjd/27r 


-100 GHz 


-12 GHz 


-6 GHz 


6 GHz 


12 GHz 


100 GHz 








(n, L) = (36, 34) state 








a«(W2)(36,34) 


-602 


-5.02 x 10 3 


-1.00 x 10 4 


1.00 x 10 4 


5.02 x 10 3 


602 


a«(wi)(36,34) 


-1.61 


-1.61 


-1.61 


-1.61 


-1.61 


-1.61 


<*t(W2) (36,34) 


601 


5.02 x 10 3 


1.00 x 10 4 


-1.00 x 10 4 


-5.02 x 10 3 


-603 


«t(wi) (36,34) 


0.204 


0.203 


0.203 


0.203 


0.203 


0.203 


Ps (W2) (36, 34) 


-0.403 


-0.401 


-0.401 


-0.401 


-0.401 


-0.400 


Pt (W 2 ) (36, 34) 


-0.945 


-0.946 


-0.946 


-0.947 


-0.947 


-0.948 








(n, L) = (34, 32) state 








a s (W2)(34,32) 


-2.06 


-2.06 


-2.06 


-2.06 


-2.06 


-2.06 


««(wi)(34,32) 


-459 


-3.81 x 10 3 


-7.62 x 10 3 


7.62 x 10 3 


3.81 x 10 3 


456 


«t(W2)(34,32) 


0.265 


0.265 


0.265 


0.265 


0.265 


0.264 


Ot(wi) (34,32) 


416 


3.47 x 10 3 


6.94 x 10 3 


-6.95 x 10 3 


-3.48 x 10 3 


-417 


Ps (Wl) (34,32) 


-1.36 


-1.36 


-1.36 


-1.36 


-1.36 


-1.36 


A(wi)(34,32) 


-0.401 


-0.401 


-0.401 


-0.401 


-0.402 


-0.402 
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=(36,34) 



Ata d 



FIG. 1: Energy level diagram indicating the two-photon tran- 
sition (n, L) = (36, 34)->(34, 32) in p 4 He + . The relative posi- 
tion of the virtual intermediate state for a detuning frequency 
Auid of the two counter-propagating laser beams are shown. 




L=29 30 31 32 33 



FIG. 2: Portion of the energy level diagram of p 4 He + . The solid lines indicate radiation-dominated metastable states, the 
wavy lines Auger-dominated short-lived states. The broken lines indicate pHe 2+ ionic states formed after Auger emission, the 
curved arrows indicate Auger transitions with minimum \AL\. Calculated values of the dipole moments for some single-photon 
radiative transitions are shown in atomic units. 




FIG. 3: Transition matrix element of the two-photon reso- 
nance (n,L) = (36, 34)— > (34, 32) as a function of the offset 
AuJd/2-K of the two laser beams from the virtual interme- 
diate state. The matrix element is averaged over the mag- 
netic quantum number M. The arrow indicates the position 
wherein the two laser beams have equal frequencies. 




FIG. 4: Scalar dipole polarizability of states (n, L) — (34, 32) and (36, 34) versus frequency offset Ao; [ j/27r of the virtual inter- 
mediate state from the real state (35,33): a s (wi) (34,32) and Qs(w2)(36,34) as defined in Eq. (J5J (a)-(b). Tensor polarizabilities 
at (wi) (34,32) and a t {w 2 ) (36,34) (c)-(d). 
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FIG. 5: Residual scalar dipole polarizability of states (n, L) = (34,32) and (36,34) versus frequency offset Aujd/2-K of the 
virtual intermediate state from the real state (35,33): /3 s (wi) (34,32) and /3 s (wi) (34,32), and /3t(w2)(36,34) and fit{w2) (36,34). The 
contributions of all states except the intermediate one (35, 33) are included. 
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FIG. 6: The ac Stark shift in the p 4 He resonance (36, 34)— ^(34, 32) estimated from the state polarizabilities a 3 and at, as a 
function of the intensity ratio I2/I1 between the two laser beams with I\ fixed at 5 x 10 4 W/cm 2 . The shift for transitions 
involving states involving the magnetic substates M = 0, ±16, and ±32 are indicated. The virtual intermediate state is tuned 
Acj,j/27r = ±12 GHz (a) and ±6 GHz (b) from the real state. 
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FIG. 7: Expected positions and intensities of the M-sublines of the two-photon transition (36, 34)— >(34, 32) in p*He , wherein 
the virtual intermediate state is offset Aujd/2n ~ 12 GHz (a)-(c) and —12 GHz (d)-(f) from the state (35,33). The intensity 
ratios I2/I1 between the two laser beams were varied between 0.6, 0.78, and 1.0 in each figure, whereas I\ was kept constant 
at ~ 5 x 10 4 W/cm 2 . 
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FIG. 8: Energy level diagram showing the hyperfine sublevels of the resonance parent, intermediate, and daughter states 
involved in the two- photon transitions of f^He 4 ^ (a). The spin orientation (S e ,Sp) of each hyperfine sublevel is indicated by 
arrows. The four strongest two-photon transitions are indicated with solid arrows. The same figure in the case of p°He ' , 
showing the eight strong lines between hyperfine levels (S e , Sh, Sp), see text. 




FIG. 9: Depopulation efficiency of the two-photon transition 
(n, L) = (36, 34)— ^(34, 32) in p 4 He + at various intensities of 
the two laser beams. 
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FIG. 10: Simulated resonance profiles of five two-photon transitions in p 4 He + excited at laser intensities p = 1 mj/cm 2 . 
The virtual intermediate state was tuned Aw d /27r ~ — 12 GHz from a real state. The positions of the four hyperfine lines 
are indicated by arrows, together with the principal and angular momentum quantum numbers [n, L) and spin orientations 
(S e , Sp) —> (S' e , Sy) of the resonance parent and daughter states. 
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Laser detuning (GHz) Laser detuning (GHz) 

FIG. 11: Simulated resonance profiles of two two-photon transitions in p 3 He + excited at laser intensities p = 1 mj/cm 2 . 
The virtual intermediate state was tuned Awd/27r ~ — 12 GHz from a real state. The positions of the eight hyperfine lines 
are indicated by arrows, together with the principal and angular momentum quantum numbers (n, L) and spin orientations 
(S e , Sh, Sp)—¥(S' e , S' h , Sip) of the resonance parent and daughter states. 
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Laser detuning (GHz) Laser intensity ratio l 2 / 1, 

FIG. 12: Simulated profiles of the resonance (n, L) = (36, 34)— >(34, 32) in Ji 4 He + , with the frequency offset Au)d/"2iT ~ —12 
GHz and intensity ratio between the two lasers I2/I1 = 0.5 (dashed lines) and 1 (solid lines) for constant I\ ~ 5 x 10 4 W/cm 2 
(a). The ac Stark shift in the simulated profiles for offsets AuJd/2-K ~ — 12 GHz (solid circles) and 12 GHz (squares), as a 
function of I2 /h ■ 



Calculation of transition probabilities and ac Stark shifts in two-photon laser 

transitions of antiprotonic helium 
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Numerical ab initio variational calculations of the transition probabilities and ac Stark shifts 
in two-photon transitions of antiprotonic helium atoms driven by two counter-propagating laser 
beams are presented. We found that sub-Doppler spectroscopy is in principle possible by exciting 
transitions of the type (n, L)—>(n — 2, L — 2) between antiprotonic states of principal and angular 
momentum quantum numbers n ~ L — 1 ~ 35, first by using highly monochromatic, nanosecond 
laser beams of intensities 10 4 — 10 5 W/cm 2 , and then by tuning the virtual intermediate state 
close (e.g., within 10-20 GHz) to the real state (n — 1, L — 1) to enhance the nonlinear transition 
probability. We expect that ac Stark shifts of a few MHz or more will become an important source 
of systematic error at fractional precisions of better than a few parts in 10 9 . These shifts can in 
principle be minimized and even canceled by selecting an optimum combination of laser intensities 
and frequencies. We simulated the resonance profiles of some two-photon transitions in the regions 
n — 30-40 of the p 4 He + and p 3 He + isotopes to find the best conditions that would allow this. 

PACS numbers: 36.10.-k, 31.15.A-, 32.70.Jz 



I. INTRODUCTION 



The transition frequencies f eX p of antiprotonic helium 
atoms [H-il (pHe + = pT + e~ + Hc 2+ ) have recently 
been measured by single-photon laser spectroscopy to a 
fractional precision of ~ 1 part in 10 8 [6j-l8j. By com- 
a ring these results with three-body QED calculations 
3|, the antiproton-to-electron mass ratio has been 
determined as 1836.152674(5) @, 0]. To further in- 
crease the experimental precision on is cxp , we have pro- 
posed future experiments [l5j of sub-Doppler two-photon 
spectroscopy of pHe + by irradiating the atom with two 
counter-propagating laser beams [16j . Dynamic (ac) 
Stark effects are expected to become one of the important 
sources of systematic error in these future experiments, 
as is the case with other high-precision laser spectroscopy 
measurements of atomic hydroge n I17H221 and antihydro- 
gen dHim, molecular hydrogen |25ll26|. helium [27H30| . 
and muonium [3lH33j |. In this paper we calculate the 
transition probability and ac Stark shift involved in these 
two-photon transitions using precise three-body wave- 
functions of pHe + . 

The pHe+ atoms can be easily synthesized by sim- 
ply allowing antiprotons to slow down [34T - I391 ] and come 
to rest in a helium target. Some of the antiprotons are 
captured [4fj| - |49j into Rydberg pHe + states with large 
principal (n ~ 38) and angular momentum [L ~ n — 1) 
quantum numbers that have microsecond-scale lifetimes 
against antiproton annihilation in the helium nucleus. 
The longevity is due to the ground-state electron in pHe + 
which protects the antiproton during collisions with other 
helium atoms [H(3-[52j . All laser spectroscop y ex periments 
Q reported so far have used pulsed lasers [53} to induce 
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FIG. 1: Energy level diagram indicating the two-photon tran- 
sition (n, L) = (36, 34)->(34, 32) in p i Ue + . The relative posi- 
tion of the virtual intermediate state for a detuning frequency 
Aula of the two counter-propagating laser beams is shown. 

single-photon transitions of antiprotons occupying these 
metastable states, to short-lived states with nanosecond- 
scale lifetimes against Auger emission of the electron (541 - 
[58|. A Rydberg pHe 2+ ion [59|-[62j then remained after 
Auger decay, which was rapidly destroyed by collisional 
Stark effects. The resulting resonance profiles of pHe + 
had Doppler widths Aujd/2tt ~ 0.3 — 1.2 GHz corre- 
sponding to the thermal motion of pHc + in the experi- 
mental target at T ~ 10 K. This broadening limited the 
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L=29 30 31 32 33 



FIG. 2: (Color online) portion of the energy level diagram of f^He . The solid lines indicate radiation-dominated metastable 
states, the wavy lines Auger-dominated short-lived states. The broken lines indicate pHe 2+ ionic states formed after Auger 
emission, the curved arrows Auger transitions with minimum |AL|. Calculated values of the dipole moments for some single- 
photon radiative transitions are shown in atomic units. 



experimental precision on v exp . 

The first-order Doppler broadening can in principle be 
reduced (Fig. [T]) [l5| by irradiating pHe + atoms with two 
counter-propagating laser beams of angular frequencies 
uji and u)2 and inducing, e.g., the two-photon transition 
(n, L) — > (n — 2,L — 2). This results in a reduction of 
Awe by a factor \(uii —U)2)/(oji +0J2)\- Among a number 
of possible two-photon transitions, a particularly strong 
signal is expected for (n, L) = (36, 34) -> (34, 32HFig. QJ, 
as this involves a large antiproton population [40 . |4l| in 
the resonance parent state (36,34). Whereas the states 
(36,34) and (35,33) are metastable with 1-^s-scale life- 
times, the resonance daughter state (34,32) is Auger- 
dominated with the lifetime r ~ 4 ns [g, O 13, HH that 
corresponds to the natural width of a spectral line of ~ 30 
MHz. 

The probability of inducing the transition can be en- 
hanced [l5l . [63l - [6g | by tuning uj\ and 0J2 so that the virtual 
intermediate state of the two-photon transition lies close 
(e.g. \AuJd/2n\ < 10-20 GHz) to the real state (35,33) 
such that, 

Wl(SI) = 47ri?ooC [-£(35,33) ~ #(34,32) (a.u.)] + Auj d , 

(1) 

W 2 (SI) = inRooC [#(36,34) - £(35,33) ( a - u -)] - Au} d . 

Here R^c = 3.289842 x 10 15 Hz denotes the Rydberg 
constant, and £(„.l) the binding energy of the p 4 He + 
state (n, L) in atomic units. At experimental conditions 
where this offset is much larger than the Doppler width 
flAcddl » Awn), the two-photon transition is expected 
to directly transfer the antiprotons populating the par- 
ent state (36,34) to the daughter state (34,32), whereas 
the population in the intermediate state (35, 33) will be 
unaffected. 



This paper is organized in the following way. Some 
details of the numerical methods are described in Sec- 
tion [H] The transition amplitude of the two-photon res- 
onance (36, 34)— >(34, 32) of |5 4 He, and the polarizabilities 
of the parent and daughter states (36, 34) and (34, 32) at 
various offsets Aw^ are estimated in Sections IIII Al and 
IIIIBI We next calculate the "background" polarizabil- 
ity due to the contributions of J5He + states other than 
the resonant intermediate state (35,33) (Section IIII CI) . 
Based on these results, the ac Stark shift and broaden- 
ing are semi-analytically estimated in Section IIII Dl We 
next discuss the hyperfine structure in two-photon tran- 
sitions of the p 4 He + and p 3 He + isotopes (Section IIII E[) 
[77M82I ] . We numerically simulate the profile of several 
two-photon resonances (Section IIII F[) before concluding 
the paper. Atomic units (a.u.) are used to evaluate the 
state polarizabilities and transition amplitudes (70l - [75| . 
whereas International System of Units (SI) are used for 
the transition frequencies, rates, and laser intensities rel- 
evant for future spectroscopy experiments. 



II. DETAILS OF THE CALCULATION 

For simplicity we take the linearly polarized laser field 
aligned along the z-axis, such that the perturbation 
Hamiltonian in a laser field of frequency u> and ampli- 
tude F has the form, 

H' = -e F(t) d, F(t) = e z F cos(iui). (2) 

Here d = XL=i ^0^0 is the electric dipole moment op- 
erator. The second order correction to the unper- 
turbed eigenenergy Eq of a pHe + state vector |0) may 



3 



then be expressed as 
1 



E^ 2 



a z /(io,M 2 )-F 2 



(3) 



where (uj) is a tensor of the dynamic dipole polariz- 
ability: 



E 

9 



(0\d*\q)(q\d3\0) | <0|rf*| g >< g |rfJ|Q> 
' E - E q - uj 



Eq — E Q 



(4) 



The energy of a pHe state vector \q) is denoted by E q , 
and the summation of q is over all states which are ac- 
cessible via single-photon transition from the resonance 
parent state |0). 

The tensor (uj) may be rewritten in terms of the 
irreducible scalar and tensor polarizability operators, 



a. 



= i[a 
3 L 



- + a Q 
a+L(2L-l) 



VV - -L 2 

3 



a (2L-l) 



(5) 



3(L+l)(2L + 3) 3(L+1) 



where the angular momentum operator is denoted by L. 
The coefficients a + , ao, and a_ are defined as follows, 



a+ = 



a 



2L 



2 _^(^-£7,)|(0L||d||?(£+l))| 

1 q 



(e - E q y - ^ 



(E -E q )\(OL\\d\\qL)\ 



2L 



1 ^ 



(Eo - E q y 



(6) 



2L 



2 _^(£fc-£,)|(0L||d||g(£-l))| 

9 



(^0 



Here a+ and a_ represent the contributions from antipro- 
ton transitions to states of normal parity which change 
the orbital angular momentum quantum number L of the 
antiproton by 1 or — 1 . The contribution ao involves tran- 
sitions to pHe + states of anomalous parity in which the 
L-value is unchanged and the ls-electron is excited to, 
e.g., the 2p state. 

For our analysis it is convenient to define " background" 
polarizabilities using the above equations, where the 
dominant contribution from the intermediate state of the 
two-photon transition is subtracted. For example, the 
corresponding scalar and tensor background polarizabil- 
ities of state (n, L) — (36, 34) can be calculated as, 



2(E -Ej) |(0L||d||z(X-l))| 
3(2L + 1) 

2(E -Ei) 



(E - EiY - oj 2 
|(0£||d|K(£-l)}| 



(7) 



3(2L + 1) (Eo-Ei) 2 



Here E and Ei denote the energies of state (36, 34) and 
the intermediate state (35,32). 



The transition matrix element hl.l-i,m of the pHe + 
two-photon transition (n, L, M )— 2, L~ 2, M ) induced 
by two linearly-polarized laser beams of total frequency, 
wi + W2 ~ 47ri? 00 c(i?o —E\), can be calculated using the 
Wigner 3j-symbols as, 



L-l 1 L 
M —M 



L-l 1 L-2 
M —M 



E 



(l(L-2)||d||g(L-l))(g(L-l)||d||0L) 



Eq — Eg 



Ll>2 



(8) 



(l(L-2)||d||«(L-l))( 9 f(L-l)||d||0L> 



Eq — E q 



wherein |1) denotes the state vector of the resonance 
daughter state of energy Ex. This is related to the two- 
photon Rabi oscillation frequency (in atomic units) of 
this laser transition via the equation, 



^2 7M (a.u.) = i \kl,l-2,m\ FxF 2 . 



(9) 



The last term in Eq. ([5]) can be neglected at small offset 
frequencies Avid- 
la order to calculate these quantities we must evaluate 
the reduced matrix elements for the dipole operator d 
and diagonalize the Hamiltonian. For this we use the 
variational exponential expansion described in Ref. Q. 
The wave function for a state with a total orbital angular 
momentum L and of a total spatial parity it = (— 1) L is 
expanded as follows, 

*£ M (R, ri )= E 3l 1 ^(R,? 1 )Gfj a (Ji,r 1 ,rs), 

h+l 2 =L 
N 

Gfj 2 (i?,r l7 r 2 ) = {c n Re[e- a " R -^-^] ( 10 ) 



+D n Im[e- Q " fl - /3 " ri - 7 " r2 ]}. 



where the complex exponents a, j3, and 7 are generated 
in a pseudorandom way, R and r! are position vectors of 
an antiproton and an electron with respect to a helium 
nucleus, and ri the distance between the antiproton and 
electron. Further details may be found in Refs. [9l. Hfjj. 

This method has been previously employed to calcu- 
late the nonrelativistic energies of_»He + with a relative 
precision of around 1 part in 10 12 [9Hllj. We here deter- 
mined the nonrelativistic values of the transition matrix 
elements using the same wavefunctions to a precision of 
~ 10 ~ 6 . This was more than adequate for our purpose of 
roughly estimating the two-photon transition probabili- 
ties and ac Stark shifts relevant to future spectroscopy 
experiments. The pHe + states of unnatural (or anoma- 
lous) parity (— 1) L+1 involve an electron in an excited 
state and therefore are not metastable (76|. They will 
not be considered here. 
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III. RESULTS 



Transition matrix element 



We first calculated the transition amplitude kl,l-2.m 
for the two-photon resonance (n,L) = (36, 34)— s>(34, 32) 

in p 4 He + at various offsets ujd from the intermediate 
(35, 33) state, and estimated the laser intensities needed 
to drive this transition. In Fig. [21 sequences of single- 
photon transitions connecting the states (36, 34) and 
(34, 32) are indicated by straight arrows, together with 
the corresponding dipole moment |(0L||d||g(L — 1)}|. The 
atomic units shown here can be converted to SI units, the 
corresponding spontaneous decay rates in s _1 obtained 
using the equation, 



7 (SI) 



e 2 a 2 H 3 |(OL||d|| g (£~l))| 2 (a.u.) 
47reo She 3 2L + 1 



(11) 



The Si-unit constants that appear in the above equation 
are, e: the elementary charge, clq: the Bohr radius, eo: 
the dielectric constant of vacuum, Ti: the reduced Planck 
constant, and c the speed of light. The angular transition 
frequency between states \q) and |0) are denoted by uj q0 . 

Two types of transitions (n, L) — > (n, L — 1) and 
(n, L)— >(n— 1, L — 1) have the largest amplitudes of ~ 1 
a.u., but the latter kind which conserve the vibrational 
quantum number v — n — L — 1 and involve fluorescence 
photons of frequency Lu q o/2ir ~ 10 15 Hz constitute the 
dominant channels of spontaneous decay. These tran- 
sitions are most favorable for laser spectroscopy. The 
transition frequencies, dipole moments, and decay rates 
of some single-photon resonances in p 4 He + and |5 3 He + of 
the type (n, L)— >(n— 1, L— 1) are shown in TableHl The 
dipole moments for the higher-lying infrared transitions 
involving states with n ~ 40 are relatively large (> 2 
a.u.), whereas for UV transitions in the n < 33 regions it 
is reduced to < 1 a.u. On the other hand, the radiative 
decay rates increase for lower-n transitions, e.g. from 
4.7 x 10 5 s- 1 for (40,36)^(39,35), to 6.6 x 10 5 s" 1 for 
(32, 31)— t-(31, 30) due to the w 3 -dependence. 

Using the single-photon dipole moments calculated 
above, we derive the two-photon transition amplitude 
k-l.l-2,m of the resonance (36, 34)— )>(34, 32) inp 4 He + for 
cases where the virtual intermediate state is offset over a 
large range between Awd/27r = —0.6 and 0.6 PHz from 
the state (35,33). The atom is excited by two linearly- 
polarized, counterpropagating laser beams. Fig.[3](a)-(c) 
show the amplitude |«£ £— 2] averaged over all ~ 70 tran- 
sitions between the magnetic substates which conserve 
the M- value, 
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FIG. 3: Transition matrix amplitude of the two- 

photon resonance (n, L) — (36, 34)— »(34, 32) as a function of 
the offset Aud/2n in the two laser frequencies uui and 0J2 from 
the virtual intermediate state. Transition probability is aver- 
aged over the magnetic quantum number M, see text. Arrow 
indicates the position where the two laser beams have equal 
frequencies. 



than the amplitude ~ 7.8 a.u. 20] for the ls-2s two- 
photon transition of atomic hydrogen excited by 243- 
nm laser light. Gigawatt-scale laser intensities would be 
needed to induce the antiprotonic transition within the 
microsecond-scale lifetime of pHe + . On the other hand, 
the transition probabilities can be strongly enhanced to 
3> 10 3 a.u. by tuning the virtual intermediate state 
within ~ 20 GHz of the real states (n,L) — (34,33), 
(35,33), or (36,33). Eq. [S] indicates that the transition 
can then be induced using nanosecond laser pulses of elec- 
tric field F ~ (1 — 2) x 10~ 6 a.u. According to the equa- 
tion, 



The lKL^-21-values are usually small, e.g. a few a.u. 
for lasers of equal frequency (u>i = u^). This is smaller 



1 f eF(a.u.) 

2 \ 47re ao 



(13) 
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TABLE I: Some single-photon transitions in p 4 He + and p 3 He + : their vibrational quantum number v, transition frequencies, 
dipole moments, and radiative rates. 



(n,L) - 


> (n',L') 


V 


Trans, freq. 


Dipole moment 


Rad. rate 








(THz) 


(a.u.) 


(10 5 s _1 ) 








p 4 He + states 






(40, 36) - 


-> (39,35) 


3 


444.8 


2.28 


4.71 


(39,35) - 


-)■ (38, 34) 


3 


501.9 


2.02 


5.44 


(38,35) - 


■> (37, 34) 


2 


566.1 


1.82 


6.36 


(37, 34) - 


(36, 33) 


2 


636.9 


1.58 


7.01 


(37,35) - 


-> (36, 34) 


1 


638.6 


1.61 


7.15 


(36, 34) - 


(35,33) 


1 


717.5 


1.38 


7.63 


(35,33) - 


■> (34, 32) 


1 


804.6 


1.16 


7.92 


(33,32) - 


-)■ (32,31) 





1012.4 


0.79 


7.42 


(32,31) - 


■> (31,30) 





1132.6 


0.62 


6.62 








o -f 

p tie states 






(39, 35) - 


-> (38, 34) 


3 


445.8 


2.30 


4.96 


(38, 34) - 


j. C37 33) 


3 


505.2 


2.03 


5.77 


(37, 34) - 


4 (36, 33) 


2 


572.0 


1.83 


6.81 


(36,33) - 


-)■ (35,32) 


2 


646.2 


1.58 


7.55 


(35,33) - 


-»■ (34, 32) 


1 


730.8 


1.37 


8.27 


(34, 32) - 


4 (33,31) 


1 


822.8 


1.15 


8.59 


(33,32) - 


4 (32,31) 





928.8 


0.95 


8.44 


(32,31) - 


4 (31,30) 





1043.1 


0.77 


7.94 



this corresponds to a peak intensity of I ~ 10 4 -10 5 
W/cm 2 which is achievable using titanium sapphire lasers 
of narrow linewidth [l6| . The second maximum at 
Wd/27r ~ —85 THz in Fig. [3] (b) corresponds to the case 
of the w\ and L02 lasers resonating with the respective 
transitions (36,34)^(35,33) and (35, 33)-^(34, 32). 



B. Polarizabilities 



SI units according to the equation, 



K(a.u.) = a(SI) x 



(14) 



The graphs follow a reciprocal Au^ 1 dependence and 
are approximately symmetric with respect to the origin. 
This simple behavior suggests that the ac Stark shift is 
primarily caused by the contribution from the interme- 
diate state (35, 33). 



We next evaluate the polarizabilities of the parent 
and daughter states of the transition. In Figs. |4] (a) 
and (c), the scalar and tensor polarizability components, 
a s (ui) (34,32), at (^l) (34,32): of state (n,L) = (34,32) are 
shown. To simplify the calculation, we initially assume 
that the atom is irradiated by a single laser field of fre- 
quency u>i (see Fig. [TJ corresponding to offsets between 
Auj d /2n = -20 and 20 GHz from the state (35,33). 
A similar plot for the polarizabilities a s (w2)(36,34) and 
Q!t(w2)(36,34) of state (36, 34) under irradiation by a laser 
field of W2 are shown in Figs. @] (b) and (d). 

As the laser frequencies are offset from Awc;/27r = 
— 100 to -6 GHz, the scalar polarizabilities decrease from 
~ —500 to ~ — 1 x 10 4 a.u., whereas the tensor polar- 
izabilities have opposite sign and increase from ~ 500 
to ~ 1 x 10 4 a.u. (Table |H]). These polarizabilities of 
1000-10000 a.u. correspond to 50-500 Hz/(W/cm 2 ) in 



C. Contribution of nonresonant states 

We next study the non-resonant or " background" con- 
tributions (Eq. [7]) to the polarizabilities from all pHe + 
states other than the intermediate state (n, L) = (35, 33). 
This will allow us to estimate how far the measured ac- 
Stark shift would deviate from the predictions of a sim- 
ple three-level model which include only the resonance 
parent, daughter, and intermediate states. Figs. [S] (a) 
and (c) show the background scalar and tensor polariz- 
abilities (3s (^2) (36,34) and ^(^2) (36,34) of state (36,34) 
when irradiated with a single laser field of frequency u>2 
at offsets between Aw(j/27r = —20 and 20 GHz. They 
remained relatively constant at <~ —0.40 and 0.95 a.u. 
respectively (Table HI)). Figs. 0(b) and (d) are the corre- 
sponding plots of As (wi) (34,32) and (3 t (wi) (34,32) for state 
(34, 32) irradiated with the wi-laser. They are similarly 




FIG. 4: Scalar dipole polarizability of states (n, L) — (36, 34) and (34, 32) versus frequency offset Ao; [ j/27r of the virtual inter- 
mediate state from the real state (35,33): ct s (wi)( 3i 32 ) and a s (w>2)(36,34) as defined in Eqs. (JSJ (a)-(b). Tensor polarizabilities 

at(Wl)(34,32) and a t{w 2 ) (36,34) (c)-(d). 



TABLE II: Scalar and tensor polarizabilities of the p*He + states (36, 34) and (34, 32) at offset frequencies of the lasers between 
Acu d /2ir = -100 and 100 GHz, see text. 

Polarizabilities (a.u.) 



ALJ d /27T 


-100 GHz 


-12 GHz 


-6 GHz 


6 GHz 


12 GHz 


100 GH 








(n, L) = (36, 34) state 








Q s(<^2)(36,34) 


-602 


-5.02 x 10 3 


-1.00 x 10 4 


1.00 x 10 4 


5.02 x 10 3 


602 


a s(^l)(36,34) 


-1.61 


-1.61 


-1.61 


-1.61 


-1.61 


-1.61 


Q t (^2) (36,34) 


601 


5.02 x 10 3 


1.00 x 10 4 


-1.00 x 10 4 


-5.02 x 10 3 


-603 


a t (Ul) (36,34) 


0.204 


0.203 


0.203 


0.203 


0.203 


0.203 


Ps (W2) (36, 34) 


-0.403 


-0.401 


-0.401 


-0.401 


-0.401 


-0.400 


/3t(^2)(36,34) 


-0.945 


-0.946 


-0.946 


-0.947 


-0.947 


-0.948 








(n, L) = (34, 32) state 








Q!s(^2)(34,32) 


-2.06 


-2.06 


-2.06 


-2.06 


-2.06 


-2.06 


Qs(wi) (34,32) 


-459 


-3.81 x 10 3 


-7.62 x 10 3 


7.62 x 10 3 


3.81 x 10 3 


456 


a*(W2)(34,32) 


0.265 


0.265 


0.265 


0.265 


0.265 


0.264 


a*(wi)(34,aa) 


416 


3.47 x 10 3 


6.94 x 10 3 


-6.95 x 10 3 


-3.48 x 10 3 


-417 


A>( W l) (34,32) 


-1.36 


-1.36 


-1.36 


-1.36 


-1.36 


-1.36 


A (Wl) (34,32) 


-0.401 


-0.401 


-0.401 


-0.401 


-0.402 


-0.402 



constant (—1.36 and —0.40 a.u.). All these background 
polarizabilities are at least three orders of magnitude 
smaller than the dominant contributions to a s and at 
arising from the intermediate state (35, 33) at small off- 
sets |Aw d /27r| < 12 GHz. 

The case of two counterpropagating laser fields of an- 
gular frequencies uj\ and ui2, and amplitudes F\ and F% 



irradiating the atom simultaneously will next be consid- 
ered. The perturbation Hamiltonian for this can be ex- 
pressed as, 

H' = -eF(f)d, 

(15) 

F(t) = e z [Fx cos(wit) + F 2 cos{w 2 t)] . 
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FIG. 5: Residual scalar dipole polarizability of states (n,L) = (36,34) and (34,32) versus frequency offset Aojd/2n of the 
virtual intermediate state from the real state (35,33): /3 a (u>i) (34,32) an d P s (uji) (34,32), and /3t(u>2) (36,34) an d /?t(ct>2)(36,34)- The 
contributions of all states except the intermediate one (35, 33) are included. 



As shown in Ref . [17j , the interference effect between the 
two laser fields can be neglected in the case of uj\ ^ u)2- 
The contribution to the ac Stark shift of state (36, 34) 
from the wi-laser [which is far off-resonance with respect 
to the upper single-photon transition (36, 34)— J>(35, 33)] 
is expressed by the scalar and tensor polarizabilitics 
a s (wi) (36,34) and a t (uJi) (36,34)- The calculated values at 
offsets between Awd/27r = —20 and 20 GHz were re- 
spectively — 1.6 and —0.2 a.u. (Table |TTJ) . The corre- 
sponding values for the daughter state a s (uj2) (34,32) an d 
c-t (W2) (34,32) were also small (—2.1 and 0.3 a.u.). 

We conclude that the two-photon spectroscopy exper- 
iment depicted in Fig. [TJ can be accurately simulated by 
a simple model involving three states interacting with 
two laser beams. Any non-resonant contribution from 
other pHe + states are at least three orders of magnitude 
smaller. 



D. ac Stark shifts of transition frequency 

The ac Stark shift in the transition frequency of the 
(36, 34)— >(34, 32) resonance can be analytically estimated 
as, 



Aw, 



,M 



2tt 



E, 



(2) 

(36,34,M) 



E, 



(2) 

(34, 32, A/) 



(16) 



pi 2 ) J p(2) 

^(36,34^) anCt ^(34,32,M) 

the parent and daughter states of magnetic substate M 
induced by the two linearly-polarized laser fields can be 



approximated as, 



E, 



(2) 

(36,34,M) 



1 

~2 
3M 2 



as(w 2 )( 36 ,34) 

1190 



2278 



O-t (^>2) (36,34) 



(17) 



E, 



(34,32,M) 



1 

~2 
3M 2 



(34,32 

1056 



2016 



a t (wi) 



(34,32) 



Fl 



In Fig. [6] (a), the ac Stark shift Aw a c,M/27r for M- 
values 0, ±16, and ±32, and two laser offsets Aw<j/27r = 
— 12 and 12 GHz obtained from the above equations are 
plotted as a function of the intensity ratio I2/I1 — f f / F\ 
between the two laser beams. Here I\ is fixed at ~ 5 x 10 4 
W/cm 2 while I 2 is scanned between (3 — 5) x 10 4 W/cm 2 . 
We find a positive shift (Aw ac > 0) at two combina- 
tions of laser offsets and intensities (Awj > 0, /1 » J2) 
and (Aw,j < 0, I\ <C h)- Conversely the shift is nega- 
tive Aa; ac < at (Aw rf < 0, h > h) and (Auj d > 0, 
h <C h)- In addition to the ac Stark shift, the tensor 
polarizabilities cause the resonance line to split depend- 
ing on the M-value of the involved states. At conditions 
of I2/I1 < 0.6 or > 1 the ac Stark shift and splitting 
can reach values of more than 5-10 MHz. At smaller off- 
sets |Awd/27r| = 6 GHz, the ac Stark shift and splitting 
become twice as large [Fig. E](b)]. 

The ac Stark shift arising from a s and the splitting due 
to at can be minimized by adjusting the laser intensities 
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FIG. 6: (Color online) The ac Stark shifts in the p 4 He + resonance (36, 34)— >(34, 32) estimated from the state polarizabilities 
a B and at, as a function of the intensity ratio I2/I1 between the two laser beams with 7i fixed at 5 x 10 4 W/cm 2 . Shifts 
for transitions involving the magnetic substates M = 0, ±16, and ±32 are indicated. Virtual intermediate state is tuned 
Aw d /27T = ±12 GHz (a) and ±6 GHz (b) from the state (35,33). 



to the values indicated by filled circles in Figs. |5] (a)-(b), 



E. Hyperflne structure 



h 
h 



s(wi) 



(34,32) 



&s(U2) (36,34) 



0.76. 



(18) 



An important point is that when the sign of the laser 
detuning is reversed, e.g., from Aud/2-K = —12 to 12 
GHz, the resulting ac Stark shift also reverses sign but 
its magnitude is the same, 



Aw ac (wi,w 2 ) = -Acj ac (wi -2Auj dl uj2 + 2Auj d ). (19) 



This means that the residual ac Stark shift can be can- 
celed by comparing the two-photon transition frequencies 
measured at laser offsets of opposite sign but the same 
absolute value, i.e., Aud and -Aw,;. 

Fig. [7] shows the ac Stark shifts and transition ampli- 
tudes \kl,l-2,m\ of all magnetic sublines between M = 
—32 and 32 of this two-photon resonance at laser offsets 
Auj d /2TT = 12 GHz (a)-(c) and -12 GHz (d)-(f). The 
shifts were calculated at three ratios of the laser inten- 
sities h/h = 0.6, 0.78, and 1.0 with I x fixed at 5 x 10 4 
W/cm 2 . The ac Stark effect causes a triangular profile 
with the M = transitions being strongest and shift- 
ing the most. In actual experiments involving pulsed 
laser beams, these shifts and splittings would smear out 
due to nanosecond-scale changes in the light field in- 
tensities and Ii{t)\ the observed intensities of the 
magnetic sublines would have a nonlinear dependence on 
\kl,l-2,m | 2 i h(t), and I^it) as simulated in Sec. MI Fl 



We next study the hyperfine lines that appear in the 
two-photon resonance profile. The hyperfine substates 
of the parent, intermediate, and daughter states (n,L), 
(n — 1,L — 1), and (n — 2, L — 2) in p 4 He are shown 
schematically in Fig. [8] (a) . Due to the dominant inter- 
action between the electron spin S e and the antiproton 
orbital angular momentum L, a pair of fine structure 
sublevels of intermediate angular momentum quantum 
number F = L ± 1/2 and a splitting 10-15 GHz arise. 
The interactions involving the antiproton spin S-p cause 
each fine structure sublevel to split by a few hundred 
MHz into pairs of hyperfine sublevels of total angular 
momentum quantum number J = F ± 1/2. In Fig. [5] 
(a), the spin orientations of the electron and antiproton 
(S e ,Sp) are indicated for the four hyperfine sublevels. 
For example the energetically highest-lying component 
consists of a spin-down electron and spin-up antiproton, 
i.e. (S e ,Sp) = (!t). 

In the J5 3 He + case [Fig.|8](b)], the electronic fine struc- 
ture sublevels oi F = L ± 1/2 are similarly split by ~ 10 
GHz. Each fine structure sublevel is then split into pairs 
of 3 He hyperfine sublevels of intermediate angular mo- 
mentum K = F ± 5 arising from the interactions involv- 
ing the nuclear spin S/j. The antiproton spin gives rise 
to eight hyperfine sublevels of total angular momentum 
J = K± i. The spin orientations of the three constituent 
particles (S e , Sh, <%•) are indicated for each substate in 
Fig. El (b). 

In Fig. |8] (a)-(b), the four and eight strongest two- 
photon transitions between the hyperfine sublevels of 
p 4 He + and p 3 He + are indicated by arrows. These tran- 
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FIG. 7: Expected positions and transition amplitudes \kl,l-2,m\ of the M-sublines of the two-photon transition 
(36, 34)— >(34, 32) in p 4 He , wherein the virtual intermediate state is offset Aujd/2ir ~ 12 GHz (a)-(c) and —12 GHz (d)- 
(f) from the state (35,33). The intensity ratios I2/I1 between the two laser beams were varied between 0.6, 0.78, and 1.0 in 
each figure, whereas ii was kept constant at ~ 5 x 10 4 W/cm 2 . 



sitions pass through the virtual intermediate state with- 
out nipping the spin of any constituent particle. Many 
other transitions are possible, but they all involve spin- 
flip and so their amplitudes are suppressed by three or- 
ders of magnitude or more. 



daughter states. The mixing between pairs of states in- 
duced by the lasers are denoted by p a b, Pbc, and p ac , the 
Auger decay rate of the daughter state by j c . The three 
detunings that appear in Eq. [201 can be calculated using 
the equations, 



F. Optical rate equations 

To simulate pHc + two-photon resonance profiles, we 
used the following nonlinear rate equations which de- 
scribe a three-level model, 

Opaa 



dt 

dpbb 
dt 

dpcc 
dt 

dpab 



= -lm(n 2 Mpab), 

= Im(0 2M( o a6 ) - lm(n 1M p bc ), 

= -IcPcc + Im(fiiMPi,c), 



(20) 



Q t - -tdabPab + i—^—(Paa - Pbb) + i—^—Pac, 

dpbc . , , . &1M , > . ^2M 

= -idbcPbc + l—^\Pbb - Pec) ~ i—^—pac, 

9 Pac ., , .QlM .^2M 

= -ld ac p ac + l—^—pab - l—^-Pbc- 

Here the density matrix p aa , pbb, and p cc represent the 
antiproton populations in the parent, intermediate, and 



d ab =E b -E a -(l+^ 



. 7a + lb 
c J 2 ' 



dbc = E c - E b - ( 1 J W! - 1 — - — , 

d ac = E c - E a - ( 1 + —J U2 - [I - —J U)l - 1 , 



(21) 



where v z denotes the velocity component of the atom 
in the direction of the W2-laser beam, E a , Eb, and E c 
the binding energy of the hypcrfinc states involved in the 
transition, and j a and 75 the radiative rates of the parent 
and intermediate states. The angular Rabi frequencies 
of single-photon transitions induced between the parent 
and intermediate states of magnetic quantum number M, 
and the intermediate and daughter states are respectively 
denoted by f^M and fliM- 

Values of VL^m/^ (in s _1 ) for the transition 
(0, L, F, K, J, M)-t(q, L—l, F' , K', J', M) in p 3 He + by 
linearly-polarized laser light of intensity l\ (in W/cm 2 ) 
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t ftJ=L-1/2 , 
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FIG. 8: (Color online) Energy level diagram showing the hyperfine sublevels of the resonance parent, intermediate, and daughter 
states involved in the two-photon transitions of yr'He (a). The spin orientation (S e , Sp) of each hyperfine sublevel is indicated 
by arrows. The four strongest two-photon transitions are indicated with solid arrows. The same figure in the case of p He , 
showing the eight strong lines between hyperfine levels (S e , Sh, Sp), see text. 



can be calculated as, 



^2M /c t \ /2Iieao 

— (O.l.) — 4 / — 

Ztt V £qC n 



(22) 



be derived using the Wigner 3j- and 6j-symbols, 
; 0LFKJM\ d \q(L-l)F'K'J'M ) 

J 1 J' \ 



(OLFKJM\d\q(L-l)F'K'J'M) (a.u.) 



where the transition matrix element (in atomic units) can 



;0L||d||g(i-l)> 



M —M 



x v /(2J+l)(2J' + l) |*' J r Jj 
XV(2^+1)(2X'+1)|^ \ 
x v /(2F+l)(2F' + l) ^ J 



(23) 
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FIG. 9: Depopulation efficiency in the two-photon transition 
(n,L) = (36, 34)— >(34, 32) of p 4 He + at various intensities of 
the two laser beams. The virtual intermediate state is off- 
set uid/2-K — —12 GHz away from state (35,33), so that the 
two laser frequencies coincide with the hyperfine component 
(S e ,Sp) = (ft)— Ktt) °f tne resonance line. 



Eq. (|23j) only provides approximate values for the tran- 
sition matrix elements since F and K are not exact 
quantum numbers of the three-body Hamiltonian. The 
results however agree with exact transition amplitudes 
within 1%, for all the pHe + transitions of the type 
(n, L, M)->(n - 1, L - 1, M) studied here. 

We numerically integrated Eq. to simulate the an- 
tiprotons depopulated by the two laser beams from the 
resonance parent state to the daughter state via a two- 
photon transition. Prior to laser irradiation, the an- 
tiprotons are assumed to uniformly populate the ~ 70 
M-sublevels of the parent state. The atoms follow a 
Maxwellian thermal distribution of temperature T ~ 10 
K. The temporal profiles of the laser pulses are assumed 
to be roughly Gaussian with pulse lengths At ~ 100 ns 

Fig. [5] (a)-(b) shows the efficiency e of the laser pulses 
depleting the population in the parent state of the two- 
photon transition (36,34)^(34,32) of p 4 He + [i.e. e = 1 
if the laser induces all the antiprotons occupying state 
(n, L) to annihilate, and e — when no such anni- 
hilations occur]. The virtual intermediate state is off- 
set LOd/liT = —12 GHz away from state (35,33), so 



that u>i + U2 coincides with the hyperfine component 
(S e , Sp) = (ft) - Ktt) °f tne resonance line. As the laser 
intensity is increased between p = and 0.2 mJ/cm 2 , 
the £-value increases quadratically as expected for a two- 
photon process. It begins to saturate at p > 1 mJ/cm 2 
corresponding to e ~ 0.3. Monte Carlo simulations in- 
dicate that the two-photon resonance signal of s ~ 0.3 
would be strong enough for clear detection against the 
background caused by spontaneously annihilating an- 
tiprotons [83| with a signal-to- noise ratio of > 5. Higher 
laser intensities would of course provide an even stronger 
signal, but power broadening effects then deteriorate the 
spectral resolution to several hundred MHz and so this 
should be avoided for high-precision spectroscopy. 

The resonance profiles of the two-photon transitions 
(40, 36) -> (38, 34) of the v = 3 cascade, (38, 35) -> 
(36,33) of v = 2, (36,34) -> (34,32) of v = 1, and 
(33,32) — ► (31,30) of v — in p 4 He + at temperature 
T ~ 10 K are shown in Figs. [TO] (a)-(d). These res- 
onances have among the largest transition amplitudes, 
and the Auger decay rates of the daughter states are large 
7a = (2.5 — 4) x 10 8 s _1 which is a necessary condition to 
obtain a strong annihilation signal [83j]. The intensities 
of the two lasers are around p ~ 1 mJ/cm 2 . The laser fre- 
quency u>i is fixed to an offset Au>d/2iv — — 12 GHz from 
the intermediate state, whereas u>2 was scanned between 
-0.9 and 0.9 GHz around the two-photon resonance de- 
fined by wi + cj2- In each simulated profile, the positions 
of the four hyperfine lines are indicated with arrows to- 
gether with the corresponding spin orientations (S e , Sp). 
The ~ 200 MHz linewidth of these profiles are primarily 
caused by the large Auger width of the daughter states, 
and the residual Doppler and power broadening. 

The resonance (40, 36)— >(38, 34) shows a two-peak 
structure [Fig. [TU] (a)] with a frequency interval of ~ 1.1 
GHz which arises from the dominant spin-orbit interac- 
tion between S e and L. Each peak is a superposition 
of two hyperfine lines with a few tens of MHz spacing 
caused by a further interaction between the antiproton 
and electron spins. The asymmetric structure of the pro- 
file of Fig. [TO] (a) is due to the fact that the 25-MHz 
spacing between the hyperfine lines (S e , Sp) = (tt)~Ktt) 
and (t4-)~~ Kt4) are small compared to the 75-MHz spac- 
ing between (J/f) - K4-t) an d (44-)~~K4-l) The spacings be- 
tween the hyperfine lines becomes gradually smaller for 
lower-lying transitions involving states of smaller n- and 
w-values, e.g., 0.8 and 0.5 GHz for (38, 35)— J>(36, 33) and 
(36, 34)— K34, 32). The hyperfine lines can no longer be 
resolved for the lowest— n transition (33, 32)— K31, 30) 
[Fig. [12] (d)]. The low transition probability (Table H 
of this resonance causes the small depopulation efficiency 
e < 0.1 seen here; laser intensities of p > 2 mJ/cm 2 would 
be needed to produce a sufficient experimental signal. 

We expect the two UV transitions (n, L) = 
(36, 34)^(34, 32) and (33, 32)->(31, 30) in p 4 He + to yield 
the highest signal-to-noise ratios in laser spectroscopy ex- 
periments. This is because the parent states (36, 34) and 
(33, 32) retain large antiproton populations for long pe- 
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FIG. 10: Simulated resonance profiles of five two-photon transitions in p 4 He + excited at laser intensities p ~ 1 mj/cm 2 . The 
virtual intermediate state was tuned Aojd/2n ~ — 12 GHz from a real state. Four profiles (a)-(d) involve a resonance daughter 
state with lifetimes r < 10 ns against Auger decay. They are simulated assuming a temperature T ~ 10 K of the atom. The 
narrow resonance (e) involves a daughter state of much longer (r ~ 1 /is) lifetime, and was simulated at T ~ 1.5 K. The 
positions of the four hyperfine lines are indicated by arrows, together with the principal and angular momentum quantum 
numbers (n, L) and spin orientations (S e ,Sp) —¥ (S' e ,Sp) of the parent and daughter states. 



riods t = 3-10 /is following pHe + formation [4(|. By 
comparison, cascade processes rapidly deplete the pop- 
ulations in higher n > 37 states within 1-2 (is, and so 
the associated two-photon spectroscopy signals contain a 
large background due to the spontaneously annihilating 
pHe+ atoms 0,53- 

Higher experimental precisions on v cxp may be 
achieved by cooling the atoms to lower temperature 
and by inducing two-photon transitions between pairs of 
pHe + states with microsecond-scale lifetimes. Fig.fTUlfe) 
shows the resonance (37, 35)— J>(35, 33) of p 4 He at tem- 
perature T = 1.5 K. Both parent and daughter states 
have lifetimes of r ~ 1 /xs, and so its natural linewidth 
~ 200 kHz is two orders of magnitude smaller than in 
the other resonances Figs. [TU] (a)-(d) studied here. It is 
unfortunately difficult to measure this transition exper- 
imentally, as the present detection method requires the 
daughter state to proceed rapidly to antiproton annihila- 
tion. Cooling the pHe + atoms to such low temperatures 



may be technically challenging. 

o 4. 

The profiles of two p He resonances which are ex- 
pected to yield the highest signal to noise ratios [4(| 

(35.33) ^(33,31) and (33,32)^(31,30) at temperature 
T ~ 10 K are shown in Figs. Qj] (a)-(b). The positions 
of the eight hyperfine lines and their spin configurations 
(S e ,Sh,Sp) are indicated by arrows. Due to the large 
number of partially overlapping lines, it may be difficult 
to determine the ^ CX p-values for p 3 He + with a similar 
level of precision as in p 4 He + . The problem would be 
especially acute in the case of (33, 32)— >(31, 30) [Fig. QT] 
(b)] which contain 8 sublines within a relatively small 
0.6- GHz interval. 

We finally use these numerical simulations to deter- 
mine the ac Stark shift under realistic experimental con- 
ditions. Fig. [T2] (a) shows the profiles of the resonance 

(36.34) ^(34,32) ofp 4 He + at temperature T~10 K and 
laser offsets Awd/27r = —12 GHz. They were calculated 
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FIG. 11: Simulated resonance profiles of two two-photon transitions in p 3 He + excited at laser intensities p ~ 1 mj/cm 2 and 
temperature T ~ 10 K. The virtual intermediate state was tuned Au>d/27r ~ —12 GHz from a real state. The positions of the 
eight hyperfine lines are indicated by arrows, together with the principal and angular momentum quantum numbers (n, L) and 
spin orientations (S e , Sh, Sp)—>(S' e , S' h , S^) of the resonance parent and daughter states. 
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FIG. 12: Simulated profiles of the resonance (n,L) = (36, 34)— >(34, 32) in p*He , with the frequency offset Alj,i/2-k ~ —12 
GHz and intensity ratio between the two lasers I2/I1 ~ 0.5 (dashed lines) and 1 (solid lines) for constant Ii ~ 5 X 10 4 W/cm 2 
(a). The atom is thermalized at T ~ 10 K. The ac Stark shift in the simulated profiles for offsets Auid/2n ~ — 12 GHz (solid 
circles) and 12 GHz (squares), as a function of I2/I1 (b). 



at two combinations of the laser intensities: Ji = 5 x 10 4 
W/cm 2 and I 2 = 2.5 x 10 4 W/cm 2 (broken lines) and 
h = h = 5 X 10 4 W/cm 2 (solid lines). As h/h is in- 
creased, the transition frequency shifts to larger values. 
In Fig. [12] (b), the ac Stark shifts Aw ac /27r determined 
from the simulated profiles of Fig. [12] (a) at laser offsets 
Aud/2ir = — 12 GHz are plotted using filled circles. It 
increases linearly from —2 MHz at I2/I1 = 0.5, to 5 MHz 
at I2/I1 = 1. A similar plot for offset Au!d/2n — 12 GHz 
is shown using filled squares. The two calculated sets of 



ac Stark shifts are of equal magnitude and opposite sign, 
the minimum occurring around I2/I1 ~ 0.65. 

IV. CONCLUSIONS 

We conclude that two-photon transitions in pHc + of 
the type (n, L)— >(n — 2, L — 2) can indeed be induced us- 
ing two counterpropagating nanosecond laser pulses of 
intensity ~ 1 mJ/cm 2 , for cases where the virtual inter- 
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mediate state is tuned within |Aw ( ;/27r| = 10-20 GHz of 
the real state (n— 1, L— 1). The spectral resolution of the 
measured resonances should increase by an order of mag- 
nitude or more compared to conventional single-photon 
spectroscopy. The ac Stark shifts at these experimental 
conditions can reach several MHz or more, but this can 
be minimized by carefully adjusting the relative intensi- 
ties of the two laser beams. Any remaining shift can be 
canceled by comparing the resonance profiles measured 
at positive and negative offsets iAw^ of the virtual in- 
termediate state from the real state. In practice, this 
can be done by, e.g., using a frequency comb (Hf to ac- 
curately control the frequencies u)% and ui2 of the coun- 
tcrpropagating laser beams. The UV two-photon transi- 
tions (36,34)->(34,32) and (33, 32)->(31, 30) in p 4 He + , 



and (35, 33)— !-(33, 31) in p 3 He + are expected to yield par- 
ticularly strong resonance signals that can be precisely 
measured. 
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